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We report on the observation of novel features in the minor hysteresis loops in a clean crystal of NbSe2 which displays a peak 
effect. The observed behavior can be explained in terms of a supercooling of the disordered vortex phase while cooling the 
superconductor in a field. Also, the extent of spatial order in a flux line lattice formed in ascending fields is different from 

(and larger than) that in the descending fields below the peak position of the peak effect; this is attributed to unequal degree 

of annealing of the state induced by a change of field in the two cases. 
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I. INTRODUCTION 

A variety of recent transport, magnetic and structural 
studies in weakly pinned superconductor 2H-NbSe2 sup- 
port the view that the peak effect (PE) phenomenon [1] 
in the critical current density (J c ) marks an amorphisa- 
tion of the flux line lattice (FLL) [2-10]. In a varying 
field experiment, for example, [6,9], J C (H) begins to in- 
crease from an initial low value at the onset field H^, 
reaches a maximum at the peak field H p and eventu- 
ally collapses below the measurable limit above the ir- 
reversibility field H irr ( < H C 2). The vortex matter is 
thought to undergo a transformation from a state with a 
high spatial order for H < to a highly disordered 
state for H > H p . This interpretation is usually ra- 
tionalized within the Larkin-Ovchinnikov (LO) collective 
pinning formalism [11], where the Larkin volume V c (« 
R 2 L c , where R c and L c are the radial and longitudi- 
nal correlation lengths, respectively) is a measure of the 
spatial extent of order in FLL and the critical current 
density is determined through the relation, /iqH J C (H) 
~ (n p < f 2 > /Vc) 1 / 2 , where n p and / are pin density 
and the elementary pinning force parameter, respectively. 
J c measurements therefore can reveal the extent of spa- 
tial order in the FLL. Magnetization hysteresis measure- 
ments are a convenient tool for estimating J C {H) and 
thus for detecting the occurrence of phase transforma- 
tions and associated changes in the vortex correlations 
across the PE regime. 

Recently, the PE phenomenon has received a great deal 
of attention due to a characteristically rich phenomenol- 
ogy that accompanies it. Prominent among them is the 
marked history dependence in the critical current. In 
the LO scenario, this translates into a history depen- 
dence of the Larkin volume V c , since quantities such as / 



and n p cannot be history dependent. The history depen- 
dence of the correlations implies strong metastability and 
is thus a hallmark of disorder in condensed matter sys- 
tems. Recent theoretical studies emphasizing the role of 
quenched random disorder in producing novel disordered 
(glassy) phases (see, for instance, Refs. 12-14) further 
illustrate the need to understand and ultimately unravel 
the complex effects of disorder readily observed in the PE 
regime. Ravikumar et al [9] have experimentally demon- 
strated via dc magnetization technique the presence of 
a highly disordered vortex state when a sample is field 
cooled (FC) in H < H p [6}. In crystals of 2H-NbSe 2 and 
CeRu2 with comparable levels of effective pinning, they 
had shown that the critical current density in the FC 
state is larger than that in the zero field cooled (ZFC) 
state for H < H p , i.e., J[ C (H) > J^ FC (H). They 
had also shown that the disordered FC state could be 
annealed to the ordered ZFC like state when the sample 
was subjected to a small change in magnetic field. In- 
dependently, a different type of history effects has been 
reported [15,16] in a wide variety of polycrystalline and 
single crystal samples of pure and doped C'eRu2 showing 
PE phenomenon. In these cases, the minor magnetiza- 
tion curves, starting at a field H on the forward branch 
of the hysteresis loop (such that < H < H p ), were 
reported to saturate without merging with the reverse 
leg of the hysteresis loop. The minor curves starting on 
the reverse branch, on the other hand, merge with the 
forward branch of the hysteresis loop. Roy et al [15] at- 
tribute the observed behavior to an anomalous nature 
of PE phenomenon in the mixed valent superconducting 
system CeRu2, in contrast with, and as distinct from, the 
conventional PE in 2H-NbSe2 [2] and most other weakly 
pinned low T c superconductors [9,10,17]. They propose 
[15,16] that this novel behavior of the minor loops reflects 
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a thermodynamic evidence for a first order transforma- 
tion to a new phase caused by the positional dependence 
of the order parameter in superconductors with large nor- 
mal state paramagnetism, like, heavy fermion supercon- 
ductors, mixed valent rare earth systems, etc. [18]. This 
is in contrast to the explanation based on metastability 
effects in the vortex matter caused by quenched disorder. 

Very recently, a similar effect with the minor hysteresis 
curves has been observed in a cuprate superconductor as 
well [19]. The microscopies in the cuprate system appears 
to bear little resemblance to that in the mixed valent sys- 
tems such as CeRu2- We propose that the above sets of 
experiments [9,15,16,19] exemplify the ubiquitous nature 
of thermomagnetic history dependence of J c in supercon- 
ductors in general and thus require an explanation that 
is independent of the microscopies relevant for supercon- 
ductivity in the diverse systems in which these effects are 
present. 

In this paper, through detailed measurements of mi- 
nor magnetisation curves on a clean single crystal of 
2H-NbSe 2 (belonging to the category of most weakly 
pinned samples of type II superconductors [2,4]), we 
present an understanding, based on the LO collective pin- 
ning description [11] for the observed path dependence 
in the critical currents and vortex correlations. We in- 
voke the notion of a disorder-assisted supercooling of a 
metastable disordered phase, that is otherwise thermo- 
dynamically stable only above H p . We further propose 
that changes in the applied field act as a driving force 
that helps to anneal the system, often partially. As a 
result, systems with different field histories are often in 
metastable states with different degrees of annealing and 
thus with different values of vortex correlations, leading 
to the different critical currents, as observed experimen- 
tally. In the present work, three types of isothermal mi- 
nor hysteresis loops were studied : (I) Decreasing field 
after cooling the sample in a field (FC-REV), (II) and 
(III) : decreasing /increasing field from a given point on 
the forward/reverse leg of the envelope hysteresis curve. 
Fig.l shows a schematic view of all these loops for the 
case where J C (H,T) is single- valued, independent of the 
field/temperature history. In all such cases, the minor 
curves merge into the envelope curve while always re- 
maining within it, without any overshoot effects. We 
show below that violations of this standard scenario pro- 
vide an understanding of the aforementioned anomalous 
behavior. 



II. EXPERIMENTAL 

We carried out dc magnetization measurements using 
a Quantum Design (QD) SQUID magnetometer (Model 
MPMS5) on a 2H-NbSe 2 single crystal with T c w 7.25 
K. The crystal was mounted with field applied parallel to 
its c-axis. The thermomagnetic history dependent mea- 
surements have been performed at 6.95 K, where a well 



recognized peak effect, manifesting as a sharp increase 
in the magnetization hysteresis, is observed at a rela- 
tively low field of about 1000 Oe. Occurrence of PE at 
such low fields, where the flux line lattice constant ao (= 
1600 A) is of the same order as the range of interaction 
( i.e., the penetration depth A a b) at 6.95K in this sys- 
tem [20], confirms that the sample has weak quenched 
disorder [7]. In the field range 1—2 kOe, the inhomo- 
geneity experienced by the sample in a 2 cm full-scan in 
a QD SQUID magnetometer is of the order of 0.1 Oe[21]. 
This value is much smaller than the threshold field Hu 
[22] required to change the sign of the induced shielding 
currents throughout the sample. Hu value is estimated 
by measuring the minor magnetisation curves [23,24] for 
fields above H p and was found to be about 10 Oe. Thus, 
the field inhomogeneity in a full scan of 2 cm does not in- 
troduce any error in the measured magnetization values, 
and we have recorded all the present data with a 2 cm 
scan length instead of half scan technique utilized earlier 
by Ravikumar et al [9,21], in their study of CeRu 2 and 
more strongly pinned sample of 2H-NbSe2. 

III. RESULTS AND DISCUSSION 

In Fig. 2, we show the magnetization hysteresis loop in 
the present 2H-NbSe2 crystal at 6.95 K, indicating the 
onset field of PE on the ascending field cycle (H+ t w 800 
Oe), the peak field H p («s 1000 Oe) and the irreversibility 
field Hi rr (« 1250 Oe). Within the LO collective pinning 
description, the Larkin volume V c begins to shrink at 
and FLL reaches nearly amorphous (though pinned) state 
at H p [7,10]. Note that the onset of PE is much sharper 
(see inset (ii)) on the ascending field cycle than on the 
descending field cycle. This immediately shows that at 
a given field value in the PE region, J c and thus V c in 
the vortex state are not the same on the ascending and 
descending field cycles in the PE regime. In what follows, 
we examine this non-uniqueness of the vortex correlations 
in greater detail. 

In Fig. 3 (a), we show the magnetization curves mea- 
sured for FC-REV case, i.e., magnetization of a field 
cooled sample measured in decreasing magnetic fields. 
The field value in which the sample was cooled was var- 
ied across the PE region. Note that in Fig. 3 (a), the FC- 
REV curve originating at H > H p merges with the re- 
verse magnetization envelope curve [23,24] in accordance 
with the critical state model. On the other hand, FC- 
REV curves, originating from a field H < H p , overshoot 
the reverse magnetization envelope curve in a clear de- 
parture from the conventional behavior shown in Fig.l. 
The difference between the highest magnetisation value 
recorded on the FC-REV curve and the notional equilib- 
rium magnetisation value could be taken as a measure of 
the critical current density in the FC state [25]. Thus, 
the minor curves in the FC-REV case produce higher J c 
values than those in the conventional descending part of 
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the envelope curve. 

This anomalous behavior can be explained by assum- 
ing a supercooling of the disordered phase and the an- 
nealing effect (i.e., increase in vortex lattice correlations 
or growth of Larkin domain volume V c ) due to a subse- 
quent field change. In the FC state, the FLL traverses 
through a pinned amorphous state as it is cooled down 
across the H p (T) line (see the phase diagram drawn in 
the inset (i) of Fig. 2). The highly disordered vortex state 
that is stable above PE curve, with a large density of de- 
fects/dislocations [2,5,7], is then effectively supercooled 
when the sample is cooled to a given T in a field less 
than H p (T). Vortex state obtained by cooling the sample 
in a field H is more disordered than that at the same field 
value on the descending branch of the envelope loop. In 
the latter case, the process of lowering of field (below H p ) 
induces partial annealing and produces a vortex state, 
which is more correlated and thus has a smaller critical 
current than that in the field cooled state. Furthermore, 
annealing induced by the field change as mentioned above 
is also clearly seen in the FC-REV magnetization values 
moving towards the reverse envelope curve (see inset of 
Fig.3(a)). 

The FC-REV curves initiated from different fields be- 
low H p form a family of curves. FC-REV curve originat- 
ing from H = Hp essentially retraces the reverse mag- 
netisation envelope curve below H p . As stated earlier, 
the vortex state at any field value below H p on the re- 
verse magnetization curve is the result of a gradual heal- 
ing of the disordered state existing at H = H p . Thus, the 
reverse magnetization envelope curve (for H < H p ), in 
principle need not, and in practice will not, be a mirror 
reflection of the forward magnetization envelope curve as 
the specific kinetics of the annealing processes are differ- 
ent. That the process of healing of FLL dislocations could 
continue down to a field well below (during the de- 
scending cycle) is a clear indicator of this difference. The 
data shows that for a given field H(< H p ), the Larkin 
domains on the reverse magnetization envelope curve are 
smaller than those on the forward magnetization enve- 
lope curve. In other words, J c values on the descending 
field cycle are larger than those on the ascending field 
cycle. This difference is likely to be due to the difference 
in the starting configurations on the ascending and de- 
scending branches. For the latter, one starts from a much 
more disordered state; thus the residual disorder after a 
comparable level of (incomplete) annealing is more than 
that in the former. 

Using the scenario described above, we now examine 
the minor magnetisation curves of the type II and III (cf. 
Fig.l) in the PE region of 2H— NbSe2 and compare and 
contrast them with anomalous behavior of minor mag- 
netisation curves in the PE region of CeRui [15,16]. In 
Fig. 3(b), we first show the minor magnetization curves 
of type II, initiating from different fields lying on the 
forward magnetization curve. The minor magnetization 
curves initiated from H > H p and H < HZ merge with 



the reverse magnetization curve within a field change of 
about 10 Oe. The threshold field Hji is thus estimated 
from these curves to be of the order of 10 Oe. However, 
for H~ x < H < Hp, the minor magnetization curves do 
not merge into the reverse envelope curve, because J C (H) 
values on the ascending field cycle are smaller than those 
on the descending field cycle, as asserted earlier. The sig- 
nificance of H~ t is that the disorder present at H — H p is 
maximally annealed at this field value on the descending 
cycle. 

In Fig. 3(c), we show the measured minor magneti- 
zation curves of the type III, i.e., by increasing the field 
from different points on the reverse magnetization enve- 
lope curve. These minor curves now overshoot the for- 
ward envelope magnetisation curve when the field is in- 
creased by about 10 Oe. The annealing due to the field 
increase (of about 10 Oe) is inadequate to produce a com- 
parable level of lattice order existing at the correspond- 
ing fields on the forward envelope loop. When the field 
is further increased, the residual disorder gets annealed. 
Thus the minor magnetisation curves eventually merge 
into the more ordered forward envelope curve. However, 
Roy et al[15] reported that in CeRu-2 the minor curves 
initiated from the reverse magnetisation envelope curve 
readily merge with the forward magnetisation curve, in 
contrast with our data. Whether this is a significant dif- 
ference, or merely a trivial one, in the comparative levels 
of annealing in the two instances, is unclear [26]. 

The critical current densities, j[ or and J r c ev on the 
ascending and descending field branches can now be es- 
timated from the maximum width of these two sets of 
minor magnetisation curves. We collate, in Fig. 4, the 
relative J C (H) values corresponding to three thermomag- 
netic histories of the sample, viz., the FC state and the 
states along the forward and the reverse legs of the en- 
velope hysteresis loop. The three sets of J C (H) values 
have been estimated [9,15,25] by taking the notional half 
width of the magnetization hysteresis at a given H to be 
proportional to the corresponding critical current density. 
The data in Fig. 4 can be summarized by the inequality, 

J[ C (H) > jr(H) > j£° r (H), 

which is consistent with an early report of transport data 
in a Nb crystal by Steingart et al [27]. In the framework 
mentioned above, this inequality corresponds to the least 
annealed state in the FC mode and the most annealed 
state on the forward curve, while the J r c ev is intermediate 
between the two. 



IV. CONCLUSION 

We propose that the PE phenomenon marks a true 
thermodynamic phase transformation between an or- 
dered solid that is stable below H^ (on forward envelope 
loop) and a fully disordered vortex state that is stable 
above H p . Further, this transformation is first order in 
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character. Thus, it is possible to supercool or superheat 
one phase into the regime of stability of the other phase. 
For the PE regime, when the free energies of the ordered 
and disordered phases are not significantly different, the 
metastability effects are expected to be prominent. More- 
over, thermal fluctuations are inadequate, at least below 
Hp, for the system to fully explore the phase space, which 
helps the metastability of the supercooled phase, aided by 
the fact that both phases have finite pinning. Substan- 
tial driving forces experienced when the field is changed 
allow the system to anneal (or fracture, as the case may 
be) towards the stable state. When the levels of an- 
nealing are incomplete and, as is often the case, unequal 
due to the previous thermomagnetic history, the system 
will typically exhibit different levels of correlations and 
thus different critical currents, as expected within the LO 
mechanism. The history dependence of critical current 
should thus be a generic occurrence in this regime for 
comparable levels of disorder. We emphasize that the 
proposed explanation of the various anomalous history 
dependent magnetization hysteresis data shown here and 
elsewhere [9,15,16,19] are in terms of an order-disorder 
transformation and disorder-aided supercooling. This ex- 
planation is , prima facie, independent of the specifics 
of microscopic considerations consistent with the ubiq- 
uitous nature of the phenomena under consideration. It 
is tempting to suggest that the so-called anomalous PE 
in C'eRu2 [15,16,18] may also find an explanation within 
the scenario described above; whether this is indeed the 
case remains to be concluded. 
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FIGURE CAPTIONS 

Fig. 1. Schematic behavior of minor magnetisa- 
tion curves initiated from field cooled (I) magneti- 
sation value (M FC ) and those from forward (II) 
and reverse (III) legs of the envelope hysteresis 
loop within the framework of critical state model, 
i.e., assuming J C (H) is uniquely defined by H. The 
field interval H// , corresponding to a threshold field 
change required to change the sign of shielding cur- 
rents throughout the sample, is also indicated. 

Fig. 2. Magnetization hysteresis of a NbSe2 crystal 
recorded at 6.95K for H || c. The inset (i) schemat- 
ically shows three different paths, viz., the zero 
field cooled (ZFC), field cooled (FC) and descend- 
ing fields from above H C 2- The PE line H p (T) and 
the upper critical field line H c2 (T) have been de- 
termined from the temperature dependent in-phase 
ac susceptibility data, as in Ref.3. Note that in the 
FC mode, the sample would cross the H p {T) line 
at different point each time, while reaching a given 
(H, T) value. Inset (ii) shows an enlarged view of 
the M-H loop in the PE region, indicating the on- 
set field H^, the peak field H p and the irreversibil- 
ity field Hi rr . 

Figs. 3 (a) to 3(c). Minor magnetisation curves in 
the given NbSe2 crystal at 6.95 K measured along 
three paths, as schematically sketched in Fig.l. 
The inset in Fig. 3 (a) shows the merger of FC-REV 
curves initiated from M^c(H) (where H = 0.4 kOc 
and 0.6 kOe) into the reverse envelope curve. Note 
that the minor curves in Fig. 3(b) do not readily 
overlap with reverse envelope curve, whereas those 
in Fig. 3(c) cut across the forward envelope curve. 

Fig. 4. Field dependence of J c for H | cat 6.95K in 
NbSe2 for three different histories as indicated. 



5 



-MEV (I) 




Forward 



Fig.la -Ravikumar et al 



8 



6 - 



o 



CD 

9 1 



4 - 



2 - 



- 



-2 - 



-4 - 



o 
o 







\ 



H 

f 

(H, 



t\ H c2 



-A 



O o 
o 



(i) 



ZFC 



6.8 7.0 7.2 
T(K) 




2 

^-0.4<T u - 



H 



o o 
o o 
o o 



O O o o 

cr f 



H 



irr 



P i H 




800 1200 
H (Oe) 



o o 



PE region 



NbSe 

H//c 



2 - 



6.95K - 



-6 







400 



800 



1200 1600 



H (Oe) 



Ravikumar et al- fig. 2 



0.8 



FC-REV 



0.4- 



0.0- 



-0.4- 




M FP (H) 



0.4 



0.8 



NbSe, 

H//c' 

6.95 K 



1.2 



H (kOe) 



Ravikumar et al - fig. 3a 



CO 
CO 

CX3 

o 




0.6 



0.4- 



0.2- 



0.0 z 



-0.2- 



-0.4- 



(b) 



Envelope loop 
Minor curves 



0.4 



NbSe. 
H//c 

T = 6.95KJ 




0.8 

H (kOe) 



1.2 



Ravikumar et al - fig. 3b 



0.6 



4r-0.2h 



Envelope loop 
Minor curves 



0.4 



0.8 

H (kOe) 



NbSe. 
H//c 

T = 6.95 K 




1.2 



Ravikumar et al - fig. 3c 




0.2 0.4 0.6 0.8 1.0 1.2 

H (kOe) 



Ravikumar et al Fig. 4 



